
clear
set more off


capture confirm file "Main_Results/"
if _rc mkdir "Main_Results/"

capture confirm file "Main_Results/Fig3/"
if _rc mkdir "Main_Results/Fig3/"

capture confirm file "Main_Results/Fig3/gphs/"
if _rc mkdir "Main_Results/Fig3/gphs/"

capture confirm file "Main_Results/Tables/"
if _rc mkdir "Main_Results/Tables/"



use "Data/Big_Pandemics.dta", clear
gen pandemic=1
drop year
ren endyear year
sort year
merge 1:1 year using "Data/UKGDP.dta"
sort year
merge 1:1 year using "Data/RealRateSince1310.dta", nogen
sort year
replace pandemic=0 if pandemic==.

gen pandemic100k=cond(dead>=100000,pandemic,0)

merge 1:1 year using "Data/sdeaths.dta" , nogen
gen war_deaths = sdeaths
replace sdeaths =0 if sdeaths==.
replace sdeaths= cond(sdeaths>20000,1,0)
sort year


keep year gdp gdppc  pandemic100k UKR_real sdeaths
sort year

replace gdppc = log(gdppc)*100			/*  take log */

tsset year
gen gdpgrowth = d.gdppc

cap drop b0* se0* b1* se1*  Horizon
cap drop h0* h1* l0* l1*
	gen Horizon = _n-1 if _n<=41
	gen b0=0 if Horizon==0
	gen se0=0 if Horizon==0
	gen h0=0 if Horizon==0
	gen h1=0 if Horizon==0
	gen l0=0 if Horizon==0
	gen l1=0 if Horizon==0
	gen D_b0=0 if Horizon==0
	gen D_se0=0 if Horizon==0
	gen D_h0=0 if Horizon==0
	gen D_h1=0 if Horizon==0
	gen D_l0=0 if Horizon==0
	gen D_l1=0 if Horizon==0
	gen pval_test = 0 if Horizon==0
	
sort  year
tsset year



/********** State space specification starts here *************/

constraint 1	[pot]L.pot			= 1		/* Definitional */	
constraint 2	[pot]L.gw			= 1		/* Definitional */
constraint 3	[gw]L.gw			= 1		/* Random walk assumption */
constraint 4	[gdppc]pot			= 1		/* Definitional */


cap drop zero								/* generate zero line for figures */
gen zero = 0


/************	 State space equations	 ******************/

sspace		(pot L.pot L.gw, state noconstant noerror)		/// /* Potential Output */
			(gw L.gw, state noconstant)				/// /* Growth of potential */
			(gdppc pot, noconstant ),				/// /* Output */
			covstate(identity) constraints(1/4) difficult /*method(kdiffuse)*/ iterate(25)
			
			qui{
			/* Generate smoothed estimates of potential and growth rate */
			predict y_star if e(sample), states smethod(smooth) equation(pot)
			
			}




gen rGDPUK = y_star  

gen time = year
eststo clear
forv h=0/40 {
gen drGDPUK`h' = f`h'.rGDPUK-l.rGDPUK
reg drGDPUK`h' pandemic100k sdeaths l(1/10).pandemic100k  l(1/10).sdeaths l(2/11).rGDPUK  if year>=1495
eststo ols`h'
	replace b0  = _b[pandemic100k]                      if Horizon==`h'
	replace se0 = _se[pandemic100k]                     if Horizon==`h'
	replace h0  = _b[pandemic100k]+2*_se[pandemic100k]  if Horizon==`h'
	replace l0  = _b[pandemic100k]-2*_se[pandemic100k]  if Horizon==`h'
	replace h1  = _b[pandemic100k]+1*_se[pandemic100k]  if Horizon==`h'
	replace l1  = _b[pandemic100k]-1*_se[pandemic100k]  if Horizon==`h'
	replace D_b0  = _b[sdeaths]                      if Horizon==`h'
	replace D_se0 = _se[sdeaths]                     if Horizon==`h'
	replace D_h0  = _b[sdeaths]+2*_se[sdeaths] 		 if Horizon==`h'
	replace D_l0  = _b[sdeaths]-2*_se[sdeaths]  	 if Horizon==`h'
	replace D_h1  = _b[sdeaths]+1*_se[sdeaths]  	 if Horizon==`h'
	replace D_l1  = _b[sdeaths]-1*_se[sdeaths]  	 if Horizon==`h'

	test pandemic100k = sdeaths
	replace pval_test = r(p) if Horizon==`h'
	estadd scalar pval_wars_pandemics = r(p): ols`h'
}
esttab , keep(pandemic*)
	cap drop Zero
	gen Zero=0
	twoway ///
	(  rarea h0 l0 Horizon if Horizon~=. , fc(blue)  fi(10) lcolor(white)  lpattern(none)  )  ///
		(  rarea h1 l1 Horizon if Horizon~=. , fc(blue)  fi(20) lcolor(white) lw(none)  lpattern(none)  )  ///
		(  rarea D_h0 D_l0 Horizon if Horizon~=. , fc(red%50)  fi(10) lcolor(white)  lpattern(none)  )  ///
		(  rarea D_h1 D_l1 Horizon if Horizon~=. , fc(red%50)  fi(20) lcolor(white) lw(none)  lpattern(none)  )  ///
		(  line  Zero    Horizon if Horizon~=. , lc(gray)  lw(med) lp(dash)                     )  ///		
		(  line  D_b0     Horizon if Horizon~=. , lc(red)  lw(thick) lp(dash_dot)  )  ///
		(  line  b0     Horizon if Horizon~=. , lc(blue)  lw(thick) lp(solid)  ///
				xlab(0(10)40, labsize(large) ) ylab(-5(5)15, angle(0) labsize(large) nogrid) xtitle("Years since pandemic event ended", size(large)) ///
				ytitle("Percent", size(large))) ///
		, title()    ///
		xsize(5) ysize(4) scale(1.1) ysc(r(-5 15))  scheme(s1color) ///
		legend(row(2) order(6 7 ) region(color(none)) size(vlarge)  label(6 "wars") label(7 "pandemics") ring(0) position(11) ) ///
		 name(UKrGDPUK_wars,replace) /// nodraw

	
graph save   Main_Results/Fig3/gphs/Fig3f.gph, replace
graph export Main_Results/Fig3/Fig3f.pdf, replace




/* Output a table with coefficients */
do Programs/make_table

esttab pandemics wars pval_test_wp  ///
	using Main_Results/Tables/Table1c.tex, page r ///
	se(2) b(2) star(* 0.1 ** 0.05 *** 0.01) ///
	coef(h0 "$ h = 0 $" h10 "$ h = 10 $" h20 "$ h = 20 $" h30 "$ h = 30 $" h40 "$ h = 40 $" ) ///
	mtitles("Pandemics" "Wars" "p-value") ///
	stats(Obs, fmt(0) l("Observations"))  

